A NOTE ON THE PRANDTL BOUNDARY LAYERS 



YAN GUO, TOAN NGUYEN 

Abstract. This note concerns a nonlinear ill-posedness of the Prandtl equation and an invalidity 
of asymptotic boundary-layer expansions of incompressible fluid flows near a solid boundary. Our 
analysis is built upon recent remarkable linear ill-posedness results established by Gerard- Varet 
and Dormy 2 , and an analysis in Guo and Tice [5] . We show that the asymptotic boundary-layer 
expansion is not valid for non-monotonic shear layer flows in Sobolev spaces. We also introduce a 
notion of weak well-posedness and prove that the nonlinear Prandtl equation is not well-posed in 
this sense near non-stationary and non-monotonic shear flows. On the other hand, we are able to 
verify that Oleinik's monotonic solutions are well-posed. 
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1. Introduction 

One of classical problems in fluid dynamics is the vanishing viscosity limit of Navier-Stokes 
solutions near a solid boundary. To describe the problem, let us consider the two-dimensional 
incompressible Navier-Stokes equations: 

4(::)+(."4+/«» i )(;:) + v/=,a(;: 

d x u v + d y v u = 0. 

Here, (cc, y) G T x R + and (u v , v u ) £ R x R are the tangential and normal components of the velocity, 
respectively, corresponding to the boundary y = 0. We impose the no-slip boundary conditions: 
(u u , v u )\ y=0 = 0. A natural question is how one relates solutions of the Navier-Stokes equations to 
those of the Euler equations (i.e., equations (ll.lj) with v = 0) with boundary condition v°\ y =o = 
in the zero viscosity limit? Formally, one may expect an asymptotic description as follows: 

(1-2) K\ (t, x, y) = K\ (t, x,y)+( <\ (t, x, y/y/v) 

where (u°, v°) solves the Euler equation and (u p , v p ) is the boundary layer correction that describes 
the transition near the boundary from zero velocity u u of the Navier-Stokes flow to the potentially 
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nonzero velocity u° of the Euler flow and thus plays a significant role in the thin layer with order 
0(y/v). We may also express the pressure p u as 

P u {t, x, y) = p°(t, x, y) + p p (t, x, —). 



(1.3) 



We then can formally plug these formal Ansatz into (II, ip and derive the boundary layer equations 
for (u p ,Vp) at the leading order in yjv. For our convenience, we denote Y = yj^fv and define 

u(t, x, Y) := u°(t, x, 0) + u p (t, x, Y), 

v(t, x, Y) := d y v°(t, x, 0)Y + v p (t, x, Y). 

The boundary layer or Prandtl equation for (it, v) then reads: 

dtu + ud x u + vdyu — dyU + d x P = 0, Y > 0, 
d x u + dyv = 0, Y > 0, 
u\ t =o = u (x,y) 
u\ Y =o = v\ Y =o = 0, 

limy^ +00 n = U(t,x), 

where U = u°(t,x,0) and P = P(t,x) are the normal velocity and pressure describing the Euler 
flow just outside the boundary layer, and satisfy the Bernoulli equation 

a t u + ua x u + d x p = o. 

This formal idea was proposed by Ludwig Prandtl [7] in 1904 to describe the fluid flows near the 
boundary. Mathematically, we are interested in the following two problems: 

• well-posedness of the Prandtl equation (jl,3D : 

• rigorous justification of the asymptotic boundary layer expansion. 

Sammartino and Caflisch [8j resolved these issues in an analytic setting where the initial data and 
the outer Euler flow are assumed to be analytic functions. Oleinik |6] established the existence 
and uniqueness of the Cauchy problem (|1.3p in a monotonic setting where the initial and boundary 
data are assumed to be monotonic in y along the boundary-layer profile. For further mathematical 
results, see the review paper |T|. In this paper, we address the above issues in a Sobolev setting. 
Our work is based on a recent result of Gerard- Varet and Dormy [2] where they established ill— 
posedness for the Cauchy problem of the linearized Prandtl equation around non-monotonic shear 
flows. 

In what follows, we shall work with the Euler flow which is constant on the boundary, that is, 
U = const. Also, by a shear flow to the Prandtl, we always mean that a special solution to (ll.3p 
has a form of (u s ,0) with u s = u s (t,Y). Thus, u s solves the heat equation: 

d t u s = d Y u s , Y > 0, 
u s \t=o = U s , 



(1.4) 



with initial shear layer U s , and with the same boundary conditions at Y = and Y = +oo as in 

CD. 

We shall work on the standard Sobolev spaces L 2 and H m , m > 0, with usual norms: 



1/2 m 

' M]l I,y '~ ( /„ „ \ u \ 2dxdY ) and W U Wh™ y ■= Yl Yl W d x d Y u h^ 

fc=0 i+j=k 
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For initial data, we will often take them to be in a weighted H m Sobolev spaces. For instance, we 
say «o £ e ~ aY H™ Y if e aY UQ S H™ Y and has a finite norm, for some a > (see, for example, [21 [3] 
where this type of weighted spaces is used for initial data). We occasionally drop the subscripts 
x, Y in H™ Y when no confusion is possible, and write H™ to refer to the weighted space e~ aY H™ Y ■ 

To state our results precisely, we introduce the following definition of well-posedness; here, we 
say that u belongs to U + X, for some functional space, to mean that u — U £ X . 

Definition 1.1 (Weak well-posedness). For a given Euler flow u°, denote U(t,x) = u°(t,x,Q). We 
say the Cauchy problem fjl .3j) is locally weak well-posed if there exist positive continuous 
functions T(-,-), C(-,-), some a > 0, and some integer m > 1 such that for any initial data 
Uq,Uq in U + e- aY H™ Y (T x R+), there are unique distributional solutions U\,U2 of (jl.3p in U + 

L°°(]0, T[; Hz. Y (T x M+)) with initial data Uj\t=o = Uq, j = 1,2, and there holds 

sup -u 2 (t)\\ H i 
(1.5) °^^ T 

< C(\\e aY [ul - U]\\ HT)Y , \\e aY [ul - U}\\ HTY )\\e aY [«J - 
m which T = T{\\e aY [ul - U)\\ H ^ Y , ||e^K " U]\\h™ y ). 

We note that when we choose u 2 = in the above definition, we obtain an estimate for solutions 
in the H^, Y space. We call such a well-posedness weak because we allow the initial data to be in 
H™ Y f° r sufficiently large m. 

Our first main result then reads 

Theorem 1.2 (No Lipschitz continuity of the flow). The Cauchy problem (11.3ft is not locally weak 
well-posed in the sense of Definition 



Our result is an improvement of a recent result obtained by D. Gerard- Varet and the second 
author [3] without additional sources in the Prandtl equation. In Section [5l we will show that in 
the monotonic framework of Oleinik (see Assumption (O) in Section [5]), the Cauchy problem (j 1 . 3 1) 
is well-posed in the sense of Definition 11.11 The key idea is to use the Crocco transformation to 
obtain certain energy estimates for d x u. We note that as shown in [2], the ill-posedness in the 
non-monotonic case is due to high-frequency in x and the lack of control on d x u in the original 
coordinates in (jl.3p . 

Finally, regarding the validity of the asymptotic boundary layer expansion, we ask whether one 
can write 

(;;:) (*,«,,) = (»° - f^>) (v) + (;■) a, A) + w (^) (m i), 

and 

P u (t, x, y) = {Vuyp^it, x, -y=), 

for shear flows u s and for some 7 > 0, where (M°(y),0) 4 is the Euler flow. Our second main result 
asserts that this is false in general, for all 7 > 0. Again, to state our result precisely, we introduce 
the following definition of validity of the asymptotic expansion in Sobolev spaces. 
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Definition 1.3 (Validity of asymptotic expansions). For a given Euler flow u° = u (y), denote 
U = u (0). We say the expansion (I1.6P is valid with a 7 > if there exist positive continuous 
functions T~(-, ■), C 7 (-, •), some a > 0, and some integers m! > m > 1 snc/i that for any initial shear 
layer U s in U + e~ aY Hy and any initial data uq,vq in e~ aY H™ y (1l xM + ), we can write (JL6]) 

mL oo ([0,T 7 ]; J ff ;c 1 iy (TxIR + )) withu s {0) = U s , (u u ,v u )\ t=0 = (u ,v ),' and p u e L°°([0, T 7 ]; L^ y (T x 
R+)), and there holds 

sup <C 7 (||e ay (C/ s -i7)|| Hm ,,||e^(«o^o)k- r ), 

0<t<T 7 x ' y ^ 

m w/nc/i T 7 = T 7 (||e ay (C/ s - t/)||#m' , ||e ay (tt , v )\\h™ y )• 
Our second main result then reads 

Theorem 1.4 (Invalidity of asymptotic expansions). The expansion (jl.6jl is no£ ua/id m i/ie sense 
of De finition \ 1.31 for any 7 > 0. 

We will prove this theorem via contradiction. We show that the expansion does not hold for 
a sequence of translated shear layers u Sn (t) = u SQ (t + s n ), s n being arbitrarily small, in which 
the initial shear layer u So (0) has a non-degenerate critical point as in [2]. Hence, if the Olenick 
monotone condition is violated, our result indicates that the remainder u u (t),v u (t) in the asymptotic 
expansion (1.6) can not be bounded in terms of the initial data in a reasonable fashion. 

Our result is different from Grenier's result [4] on invalidity of asymptotic expansions. He allows 
the initial perturbation data to be arbitrarily small of size v n to a fixed unstable Euler shear flow 
(could even be monotone!), and shows that in a very short time of size y^log^/V), the solution u 
grows rapidly to 0(v l / A ) in L°°. In contrast, we have to work with a family of non-monotone shear 
flow profiles, and we do not know how badly the solutions grow. On the other hand, we show that 
the expansion is invalid in order 0(u^), for any 7 > 0. In addition, the blow-up norm in [4] is H 1 
in the original variable y, whereas our result concerns the (weaker) norm in the stretched variable 
Y = y/V", upon noting that \\u(y/^/u)\\ L 2 {R+) = ^ 1/4 ||n(y)|| L 2, (IR+) . 

We remark that both our ill-posedness and no expansion theorems are in terms of a weighted H 1 
norm, which does not even control L°° norm in the two-dimensional domain. However, our proofs 
fail for a weaker space than H 1 (e.g., L 2 ), because we need the local compactness of H 1 to pass to 
various limits as v — > 0. 



2. Linear ill-posedness 

In this section, we recall the previous linear ill-posedness results obtained by Gerard- Varet and 
Dormy [2] that will be used to prove our nonlinear illposedness. For notational simplicity, we define 
the linearized Prandtl operator C s around a shear flow u s : 

C s u := —d Y u + u s d x u + vdyu s , v = — d x udy' . 

Jo 

With our notation, the nonlinear Prandtl equation (|1.3j) in the perturbation variable u := u — u s 
then reads (dropping the titles) : 

, . f d t u + C s u = -ud x u-vdYU, Y > 0, 

\ u\ t=0 = n , 

with zero boundary conditions at Y = and Y = 00. 
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Removing the nonlinear term in (|2.ip . we call the resulting equation as the linearized Prandtl 
equation around the shear flow u s : 

(2.2) d t u + C s u = 0, u\ t=0 = uq. 

Denote by T(s,t) the linearized solution operator, that is, 

T(s, t)uo := u(t) 

where u(t) is the solution to the linearized equation with u\t= s = uo- The following ill-posedness 
result is for the linearized equation f|2 . 2[) . 



Theorem 2.1 ([2]). There exists an initial shear layer U s to Q1.4)l which has a non- degenerate 
critical point such that for all Eq > and all m > 0, there holds 

(2.3) sup \\T(a,t)\\ C f H mja) = +oo, 

0<s<t<e 

where \\ ■ \\c(H m ,L 2 ) denotes the standard operator norm in the functional space C(H™, L 2 ) consisting 
of linear operators from the weighted space space. 

Sketch of proof. In fact, the instability estimate ()2.3p stated in [2] was from the weighted space H™ 
to another weighted space H™ . From their construction, (|2.3j) remains true when the targeting 
space is not weighted. We thus sketch their proof where it applies to the usual L 2 space as 
stated. We recall that the main ingredient in the proof is their construction of approximate growing 
solutions u £ to (|2.2p such that for all small e, u e solves 

8 t u £ + C s u £ = e M r £ , 

for arbitrary large M, where u £ and r £ satisfy: 

ce <Wv^ < \\ u £ (t)\\ L 2 < Ce e ° t/V£ , \\e aY r £ {t)\\ H m < Ce~ m e e ° t/V£ , 

for all t in [0, T], m > 0, and for some 9q, c, C > 0. 

The proof is then by contradiction. That is, we assume that \\T(s, Oll/Zf//™,!/ 2 ) is bounded for all 
< s < t < eq, for some > and some m > 0. We then introduce u(t) := T(0,t)u £ (0), and 
v = u — u £ , where u £ is the growing solution defined above. The function v then satisfies 

(2.4) d t v + C s v = -e M r £ , v\ t=0 = 0, 

and thus obeys the standard Duhamel representation 

r-t 

M 



v {t) = -e M / T(s,t)r £ (s)ds. 
Jo 

Thus, thanks to the bound on the T(s,t) and the remainder r £ , we get that 

\\v{t)\\ L 2 < Ce M f \\e y r £ (s)\\ H ™{s)ds < C e M+ ^~ m e^ . 
Jo 

Also, from the definition of u(t), we have 

||n(t)|| L2 = ||T(0,tK(0)|| i2 < C\\e aY u £ (0)\\ H m < C e~ m . 
Combining these estimates together with the lower bound on u £ (t), we deduce 

Ce- m > \\u(t)\\ L 2 > \\u £ (t)\\ L 2 - \\v(t)\\ L2 > (c- Ce M+l 2~ m ) e^ . 



This then yields a contradiction for small enough e, M large, and t = K\ \i\.e\s/e with a sufficiently 
large K. The theorem is therefore proved. □ 
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Next, we also recall the following uniqueness result for the linearized equation. 
Proposition 2.2. (^5\) Let u s = u s (t,y) be a smooth shear flow satisfying 

sup ( sup | u s | + / y\d y u s \ 2 dy\ < +00. 
t>0 ^ y>0 Jo > 

Letu G L°°(]0, T[; L 2 (T x K+)) with d y u G L 2 (0, T x T x M + ) 6e a solution to the linearized equation 
of (12. 2D around the shear flow, with u\t=o = 0. Then, u = 0. 

Proof. For sake of completeness, we recall here the proof in [3]. Let w G L°°(]0, T[; L 2 (T x 
and Sytf; G L 2 (0,T x T x R + ) be a solution to the linearized equation of (|2.2p around the shear 
flow, with ui|f=o = 0. Let us define w k (t, y), k G Z, the Fourier transform of x, y) in x variable. 
We observe that for each k, w k solves 

{d t w k + iku s w k - ikdyU s w k (y')dy' - dyW k = 
w k (t,0) = 
w k (o,y) = 0. 

Taking the standard inner product of the equation (12. 5p against the complex conjugate of w k 
and using the standard Cauchy-Schwarz inequality to the term f v w k dy' ', we obtain 

1 

2 



ydt\\w k \\ 2 L2(R+) + \\d y w k \\ 2 L 2 iR+) < \k\ \u s \\w k \ 2 dy + \k\ \d y u s \y l/2 \w k \ \\w k \\ L -2^dy 

< \k\( sup|u s | + / y\dyU s \ 2 dy) \\w k \ 
K t,y Jo ' 



Il2(R + )- 



■y 

Applying the Gronwall lemma into the last inequality yields 

\\Mt)hHR + ) <Ce c ^ t \\w k (0)\\mR + ), 

for some constant C. Thus, w k (t) = for each k G Z since w k (0) = 0. That is, w = 0, and the 
proposition is proved. □ 



3. NO ASYMPTOTIC BOUNDARY LAYER EXPANSIONS 

In this section, we will disprove the nonlinear asymptotic boundary-layer expansion. Our proof 
is by contradiction and based on the linear ill-posedness result, Theorem 12.11 Indeed, for some 
7 > 0, let us assume that expansion (|3. 1[) is valid with 7 > in the Sobolev spaces in the sense of 
Definition II .31 for any initial shear layer U s and initial data uo,vq G e~ aY H m (T x That is, we 
can write 

p.i) (*:) (i,«, V ) = (» -«v°) ( 9 ) + (;•) («,-^ )+ w (J ff ) (m,^), 7 > 0. 

and 

p^(t,x,y) = (^) 7 ^(t,x,^), 7 >0, 

where u^y" G L°°([0, T 7 ]; F 1 (TxM + )) and ^ G L°°([0, T 7 ]; L 2 (TxM + )), for some T 7 = T 7 (||e ay (L^- 
^H^m', ||e ay ('Uo,'Uo)||H™y)- Furthermore, for some C 7 = C y (\\e aY (U s -U)\\ H p> , \\e aY (u , v )\\ H ™ Y ), 
there holds 

sup <C 7 (||e ay (C/ s -[/)||w,||e^(fio^o)||^). 

0<t<T 7 y 
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We then let (u,v) and p be the weak limits of (u u ,v v ) and p v in L°°([0, T 7 ]; H l (J x R + )) and in 
L°°([0, T 7 ]; L 2 (T x R+)), respectively, as v — > 0. Note that it is clear that T 7 ,C 7 are independent 
of the small parameter v. 



Hence, plugging these expansions into (II. lh . we obtain 

| y =o + u s )d x u v + v v {y/vd y 



d t u u + («° - u°\ y=0 + u^m" + v v {s/Dd y u Q + 9 y u s ) + 9^ - d\u v 



and 

ud t v v + i/(u° - u\ =0 +u s + (V^yu u )d x v u + (v^) 7+ Vdy^ + dyp v = v 2 dlv v + vd\v v 

We take z/ — > in these expressions. Since (u u ,v u )(t) converges to (u, v) weakly in H , the nonlinear 
terms (u u d x + v u dy)u u and {u u d x + v u dy)v u have their weak limits in L , and thus disappear in the 
limiting equations due to the factor of (y'l/) 7 . Similar treatments hold for the linear terms. Note 
that (u° — u°\ y =o)(y) = yd y u° = ^Yd y u° also vanishes in the limit. We thus obtain the following 
equations for the limits in the sense of distribution: 

d t u + u s d x u + v&yUs - dyu + d x p = 0, 

dyp = 

and the divergence- free condition for (u, v). Prom the second equation, p = p(t, x). Setting Y = +oo 
in (|3.2p and noting that (u, v) belong to the H 1 Sobolev space and u s has a finite limit asF-y +oo, 
we must get d x p = in the distributional sense. That is, the next order in the asymptotic expansion 
solves the linearized Prandtl equation: 

(3.3) d t u + C s u = 0, u\ t =0 = uq, 

with zero boundary conditions at Y = and Y = +oo, for arbitrary shear flow u s = u s (t,Y). 

Now, let u SQ be the shear flow in Theorem 12.11 such that (|2,3p holds. Thus, we have that for a 
fixed Eq > 0, m > 0, and any large n, there are s n ,t n with < s n < t n < Eq and a sequence of Uq 
such that 

(3.4) \\e aY u%\\ Hm+ i = 1 and \\u n L {t n )\\ L 2 > 2n 

with it2(i) being the solution to the linearized equation ()2.2[) around u SQ with u2(s n ) = lift. 

For such a fixed shear flow u So , and fixed n, s n as in (|3.4p . we consider the expansion (|3.ip for 
^s„ = ^s (^ + s n) an d initial data (iiq", v^ 71 ) defined as 

(3.5) (u u > n ,vZ n ):=(u%,v%), with v% := - T d x u%dy' . 

Jo 

We note that since Uq is normalized, (uq U , Vq' 11 ) belongs to e~ aY H m with a finite norm of size 
independent of n. We let T 7 ,C 7 be the two continuous functions and (u u ' n ,v u,n ) be the corre- 
sponding solution in the expansion in L°°([0, T 7 ]; i7 1 (T x R+)) whose existence is guaranteed by 
the contradiction assumption, the Definition 11.31 Furthermore, there holds 

(3.6) sup ||(fi"' B (t),^»(t))||fli <C 7 (||e^K ( Sn )-C/)||^,||e« y K 5 <)||^), 

0<t<T 7 y 

in which T 7 = T 7 (||e QY (n S0 (s n ) - l/JH^w , \\e aY «, We note that thanks to ([331) and 

the fact that u SQ solves the heat equation (jl.4p . the norms of the translated shear layer u So (s n ) — U 



s 
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and the initial data are independent of n and u, and thus so are T 7 and CL. In addition, since £o 
was arbitrarily small so that (I3.4p holds, we can assume that Eq < T 7 . 

Next, let (u n ,v n ) be the limiting solutions of (u u,n , v v,n ) when v — > 0. As shown above, we then 
obtain the linearized Prandtl equation for u n with initial data u^: 

d t u n + C Sn u n = 0, u n \ t=0 = <. 
Thus, if we define u n (t) := u n {t — s n ), the above equation reads 

d t u n + C So u n = 0, u n \ t=Sn = v%, 

which, by uniqueness of the linear flow, yields u n = on [s n ,£o] an d 

\\u n {t n -s n )\\ L 2 = \\u n {t n )\\ L 2 > 2n. 

This implies that for small u, \\u n ' u (t n — s n )\\^2 > n, which contradicts with the uniform bound 
(13. 6p as n is arbitrarily large and C 7 is independent of n. The proof of Theorem 11.41 is complete. 



4. Nonlinear ill-posedness 

Again, using the previous linear results, Theorem 12. 1\ we can prove Theorem [L2] for the nonlinear 
equation (|1 .3|> . We proceed by contradiction. That is, we assume that the Cauchy problem (j 1 . 3 1) 
is (H m ,H l ) locally well-posed for some m > in the sense of Definition 11.11 Let u So be the fixed 
shear flow in Theorem 12.11 such that (|2.3p holds. By definition, (|2.3p yields that for fixed Eq > 
and any large n, there are s n ,t n with < s n < t n < so and a sequence of Uq such that 

(4.1) \\e aY u%\\ Hm = l and \\u n L {t n )\\ L 2 > n 

with n"(t) being the solution to the linearized equation (I2.2p around u SQ with u2(s n ) = Uq. We 
now fix n large. 

Next, define Vq n := u so (s n ) + 5uq , with 5 a small parameter less than 5$. Let v s,n be the solution 
to the nonlinear equation (|1.3p with v S ' n \t=o = Vq U . By the well-posedness applied to two solutions 
v s,n and the shear flow u Sn (t) := u so (t + s n ), there are continuous functions C(-, -),T(-, •) given in 
the definition such that 

ess sup \\v 5 > n (t) -u S0 (t + s n )\\ H i < C5\\e aY v%\\ H m = C 5, 
te[o,T] 

in which T = T(\\e aY (v S ' n -U)\\ Hm , \\e aY (u So (s n )-U)\\ H m) and C = C(\\e aY (v S > n -U)\\ H ™ , \\e aY (u So (s n )- 
U)\\H m )- Thanks to (14. ip and the fact that u so solves the heat equation (11. 4p . v^ 71 — U and u S0 (s n )—U 
have their norms in e~ aY H m bounded uniformly in n and 5. Thus, the functions T(-, •) and C(-, •) 
can be taken independently of n and 5. In what follows, we use T, C for T(-, •), C(-, •). 

In other words, the sequence u s,n := ^(v s,n — u Sn ) is bounded in L°°(0, T; i? 1 (T x R+)) uniformly 
with respect to 5, and moreover it solves 

(4.2) d t u 5 ' n + £ Sn u 5 > n = 5N(u 5 > n ), u 5 > n (0) = v$, 

noting that £ Sn is the operator linearized around the shear profile u Sn and is the nonlinear term: 
N(u s,n ) := —u S ' n d x u S ' n — v S ' n dYU S ' n . From the uniform bound on u S ' n , we deduce that, up to a 
subsequence, 

u 5 ' n u n L°°(0,T;H l {J x R+)) weak* as 5 0. 
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We shall show that u n solves the linearized equation ()2.2p in the sense of distribution. To see this, 
we only need to check with the nonlinear term. First, on any compact set K of R + , we obtain by 
applying the standard Cauchy inequality and using the divergence-free condition: 

r Y , r \ 1/2 



</ |0 x u*' n |dY' < C ^ 1/2 ( / |c\.u 5 ' n | 2 dF 
Jo ^ Jr+ 

and 

/ \u S ' n v S ' n \dYdx<C K [ I ' a...o„.\^ 



u 5 ' n \ / i^' n i 2 dy] (iy^x 



< 



Ck( f [ \u 5 ' n \ 2 dYdx) 1/2 ( [ f \d x u 5 ' n \ 2 dYdx 
WtJk ' ^ Jt Jr+ 



1/2 



< Ck\\u \\ h i, 

for some constant Ck depending on K. Now, from the divergence-free condition, we can rewrite 

N(u s ' n ) as 

N(u s ' n ) = -d x (u 5 > n ) 2 - dY(u S ' n v S ' n ) 
we have, for any smooth function <f) that is compactly supported in K, 



[ N(u 5 > n )cpdxdy 



< C K ^5 / (|n 5 ' n | 2 + \u s > n v s > n \ )dxdY 

< C K ^5\\u 5 > n « 2 



\m — T u ' 

as 8 — > 0, thanks to the uniform bound on u S ' n in H 1 . Here, Ck^ is some constant that depends 
on K and VF 1 ' 00 norm of (j). Thus, the nonlinearity 5N(u S ' n ) converges to zero in the above sense 
of distribution. This shows that by taking the limits of equation (I4.2p . u n solves 

d t u n +C Sn u n = 0, u n \ t=0 = u%. 



By shifting the time t to t — s n , re-labeling u n {t) := u n (t — s n ), and noting that by definition 
C Sn (t) = C so (t + s n ), one has 

d t u n + C so u n = 0, u n \ t=Sn =u%, 

that is, u n solves the linearized equation (12.2p around the shear flow u S(y By uniqueness of the 
linear flow (recalled in Proposition I2.2p . u n = on [s n ,T]. This therefore leads to a contradiction 
due to (|4.ip and the fact that the bound for u S ' n yields a uniform bound for u n and thus for u n : 

n<\\u n L (t n )\\ L 2 = \\u n (t n )\\ L 2< sup \\u n (t)\\ m < C, 

te[s n ,T] 

for arbitrarily large n, upon recalling that C is independent of n. This completes the proof of 
Theorem 11.21 



5. WELL-POSEDNESS OF THE OLEINIK'S SOLUTIONS 

In this section, we check that the Oleinik solutions to the Prandtl equation (j 1 . 3 [) are well-posed in 
the sense of Definition 11.11 Here, since now we only deal with the Prandtl equation, we shall write 
(x,y) to refer (x,Y) in (ll.3p . and use both d and subscripts whenever it is convenient to denote 
corresponding derivatives. To fit into the monotonic framework studied by Oleinik, we make the 
following assumption on the initial data and outer Euler flow: 
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(O) Assume that U(t,x) is a smooth positive function and d x U,dtU/U are bounded; the initial 
data uo(x, y) is an increasing function in y with uq(x, 0) = and uo(x, y) — > U (0, x) as y — > oo, and 
furthermore, for some positive constants 9q,Cq, 

(5.D «„ < Uln d ";° {x i , < c„. 

[7(0, x) - u {x,y) 

We also assume that all functions d y UQ, d x uo, d x d y uo are bounded, and so are the ratios dyUo/d y uo 
and dyUQdyUo/dyUo. 

We now apply the Crocco change of variables: 

. u(t,x,y) 
{t, x, y) (->■ (t, x, r)), with 77 := — — — — , 

U (t, x) 



and the Crocco unknown function: 

w(t, x, rj) 



d y u(t,x,y) 



U(t,x) ' 

The Prandtl equation (j 1 . 3|) then yields 

dtw + r]Ud x w — Ad v w — Bw = w 2 d^w, < 77 < 1, xeT 
(5.2) ^ (^9^ + ^£7 + <%£//£/) 1,^0 = 0, 

W|r,=l = 0, 

with initial conditions: w\t=o = u>o = d y uo/U. Here, 

A := {r, 2 - l)d x U +( v -l)^, B := - V d x U - 

To see how the boundary conditions are imposed, one notes that r] = and r/ = 1 correspond to 
the values at y = and y = +00, respectively. At y = +00, it is clear that w = d y u = since u 
approaches to U(t, x) as y — > +00, while by using the imposed conditions on u and v at y = 0, we 
obtain from the equation (II. 3p that 

= dlu - d x P = d y w + d x U + dtU/U = wd v w + d x U + d t U/U. 



Theorem 5.1. ([6\) Assume (O). Then there exists a T > which depends continuously on the 
initial data such that the problems (15.2|) and (I1.3P /lave a unique solution w and u on their respective 
domains, and there hold 

(5.3) 0i(l-77) < w(t,x,ri) < 9 2 (l-v), \d x w(t,x,r])\ , \d t w(t,x,r])\ < 2 (l-r/) 
/or aZZ (t,x,r/) G [0, T] x T x (0,1), and 

[ ' 1 " U(t,x)-u(t,x,y) - 2 ' " C/(t,x) " 

/or a// (t,x,y) G [0, T] x T x IR +; /or some positive constants O\,02. In addition, weak derivatives 
dtu,d x u,d y d x u,d y u,dyU are bounded functions in [0, T] x T x E + . 

Proof. In fact, the authors in [6l Section 4.1, Chapter 4] established the theorem in the case 
x £ [0, X] with zero boundary conditions at x = 0. Their analysis is based on the line method to 
discretize the t and x variables and to solve a set of second order differential equations in variable r\. 
It is straightforward to check that these lines of analysis work as well in the periodic case x £T with 
minor changes in the choice of boundary conditions. We thus omit to repeat the proof here. □ 



A NOTE ON THE PRANDTL BOUNDARY LAYERS 



11 



Using the estimates in Theorem 15 .1\ we are able to prove that 

Theorem 5.2. The Cauchy problem (jl,3p under the assumption (O) is well-posed in the sense of 
Definition \l.l\ with some constant a and some continuous functions T(-, •) and C(-, •) appeared in 
the stability estimate (|1.5|) that depend on Oq,Cq in our assumption (O). 



In the proof, we need the following lemma. 

Lemma 5.3. Under the same assumptions as in Theorem \5.1\ we obtain 
(5.5) I(t) < C7(0), < t < T, 

with 



J(t) :- 



\W\X ~ W 2x ? \WI - W 2 \ 2 



(t, x, n)dxdrj 



VO < p < 3, 



Tx[0,l] L (1-^ ' (1-VY 
for arbitrary two solutions w\,w 2 to (|5,2p . 

Proof of Lemma \5.3l We consider w\,w 2 being solutions to (|5.2p . We first note that I(t) is well- 
defined for (3 < 3 by the bounds in Theorem 15.11 that \wj\ < C(l — rj) and \wj x \ < C(l — r/). Let us 
introduce (j) = w± — w 2 . Then, (j) solves 

+ rjU (fix - A(j) v - Bcfr- (wt + w 2 )d^w 2 (t) = wfd 2 cp, < n < l,x £T 
{wi<f) v + w 2v <P)\ v =o = 0, 
(f>\v=i = 0, 

for A, B being defined as in (|5.2p . In particular, we have \A\ < C(l — rj) and \B\ < C. Multiplying 
the equation by e~ kv (j)/(l — rj) 13 and integrating it over T x (0, 1), we easily obtain 

1 d f e~ k ^\(p\ 2 



2* Ax (o,i) (i-vY 



dxdn 



Tx(0,l) 



nUtfix-Acfirj- 



(wi+W2)d v W2<j)-W 1 (f)r m 



(l-r,)P 



dxdrj. 



We treat each term on the right-hand side. Using the bounds on A, B and on Wj,d 2 Wj, it is easy 
to see that 



Tx(0,l) 



rjU<j) x - A(j) v - Bcp + (u>i + w 2 )d v W2(t) 



dxdn 



< e 



v | dxdrj + C £ 



'Tx(o.i) (l-f?)' 3 Y " f " "' ! " fc .Arx(o,i) (I-??)' 
for arbitrary small e. For the last term, integration by parts yields 



dxdrj, 



e kr, w 2 



Tx(0,l) (1 -VY 



\§ m 4>dxdrj 



e kr/ w 2 



K\(f) v \ 2 dxdrj 



/Tx(o,i) (! - VY 
f / e~ kr ^w 2 \ 

7tx(o,i) \{i-vr J 

Again, by integration by parts, we have 



-kr, w 2 



Tx{r,=0} (1 -Vf V 



cp v (pdx. 



Tx(0,l 



e kri w 2 



) dv Hi-vY 



-g ) 4> v cj)dxdrj 



1 



Tx(0,l) 



A(l-r/> 



Y)\ct)\ 2 dxdrj + 



0,, 



Tx{?j=0} 



e kv w 2 



(i- v y 



Thanks to the bounds \wj\ < C(l — rj), we have 



d 2 j e : kVw h <c- 



-ki) 



r '\{l-rjY) ~ (l-rjf 
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Collecting all boundary terms, we need to estimate 



Tx{??=0} 



k 6 kVw l U|2 | q ( w l \ c -knu\2 e kT ' w l 

k (l-n¥ m +dr >Ul- v )PJ e W (1 



(1-7?)^ 



dx. 



Note that at r, = 0, w\ 7^ and ifi^ = —w 2r) 4>. Thus, by taking sufficiently large in the above 
expression, we can bound it by 

k f 

w 2 \(f>\ 2 dx. 

Tx{ij=0} 



Combining the above estimates and choosing e sufficiently small, with noting that \ A\ < C(l — 77) < 
Cwi, we thus obtain 



(5.6) 



d_ 

dt 



e- kl i\(j)\ 2 



f e- kr >wl . l2 
„dxdr]+ / — rn\<f>r,\ dxdr] 

w 7tx(o,i) (i -vr 



+ I w\\(j)\ 2 dx < C I e V ^l H dxdr]. 

hx{ v =0} 7tx(0,1) v 1 ~ Vr 



To obtain estimates for (fr x , we take x-derivative of the equation for <j> and integrate the resulting 
equation over T x (0, 1) against e~ kr, (fi x /(l — 7]) 13 . We arrive at 
(5.7) 



1 d 



-kr)\ 



2d* Ax (0,1) (l- V )P 



^-rdxdrj 



V u 4>xx + vU x 4>x ~ A x (j)„ - A<j> x „ - B x 



Tx(o,l) 



{(wi + w 2 )d 2 w 2 ) x (j) - (u>i + w 2 )d 2 w 2 (p x ~ wl<t>xrm ~ 2 ^i^ix<A w 



(1-^ 



dxdr]. 



Similarly as in deriving the estimate (|5.6p . integration by parts and the bounds on A,B,Wj easily 
yields 



/ 



Tx(0,l) 



'/'/ 



dxdr] 



(5.8) 



< 



krj 2 



e "W 



1 

2 Ax(o,l) ^ l ~vY 
e~ k ^ x 



Or, 



Tx(0,l) 



(I-77) 



^\(f) X71 \ 2 dxdr] + C 
<f> xr) <j) x dxdr) + 



Tx(0 



\(j) x \ 2 + wj\(t> v \ 2 

i) (i - vf ' 



~ kr >dxdr] 



Tx{r;=0} 



w 1 (j) xri + 2wiw lx (, 



(j) x dx. 



Here, we note that there is a crucial factor of w 2 in front of the term |<^| 2 thanks to the bounds: 
Wj ~ (1 — rj) and \wj x \ < C(l — r,). Again, applying integration by parts to the third term on the 
right-hand side yields 



On 



Tx(0,l) 



e kv w 2 



(1 - n) 
1 

~ 2 



4> xr] 4) x dxdr] 



Tx(o,l) 



3 



e fc? % 2 



(1-7,)/ 



6 X I dxdr, + 



0„ 



Tx{?j=0} 



e ^w 2 



(1-7,)/ 



b x \ 2 dx, 
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where the last boundary term is clearly bounded by 



k 



Tx{?j=0} 



w\\4) x \ 2 dx. 



We now estimate the boundary term in f)5.8[) . We recall that at the boundary rj = 0, we have 
w\(p v = —w 2rj (j). Thus, 

wj(j) xv = Wi(-W 2v (f)x - W 2xV <P ~ Wi x (j) v ) = -Wi(w 2rl 4>x + W 2 xr l <t>) + Wi x W 2v (f>. 

That is, the normal derivative cj)^ on the boundary can always be eliminated to yield 



Tx{r?=0} L 



W 1 (f)xr 1 + 2w 1 Wi x C 



b x dx<C (\<j)\ 2 + \<j)x\ 2 )dx. 

JTx{ri=0} 



The remaining terms on the right-hand side of (|5.7p are again easily bounded by 



Tx(0,l) 



(1-7^ 



-dxdrj. 



Putting these estimates into (|5.7p . we have obtained 
(5.9) 

d f e~H 



dt 



Tx(0,l) (! ~ r ?) /3 



■dxdrj < C 



2 + |0x| 2 + hl| 2 K| 2 



Tx(0,l) 



dxdri + C / \(f)\ 2 dx. 

Tx{??=0} 



Adding together this inequality with a large constant M times the inequality (|5.6p . we can get rid 
of the boundary term and the term involving \4> v \ 2 on the right-hand side of (|5.9p and thus obtain 



(5.10) 



d_ 

dt 



M 



+ 



Tx(0,l) 



, -dxdr] < C(M) f e ~kv M \f +1^1 dxd 

(i-^r Ax (o,i) 0--vr 



The claimed estimate (|5.5p thus immediately follows from (|5.10p by the standard Gronwall inequal- 
ity, and this completes the proof of Lemma 15.31 □ 



We are now ready to give 

Proof of Theorem \5.2l We only need to check the stability estimate (jl.5p . Let U(t,x) be a fixed 
Euler flow, and take uq\{x, y) and uq 2 {x, y) be arbitrary smooth functions satisfying the assumption 
(O). Let u\,u 2 be solutions to (I1.3P and w\,w 2 the corresponding solutions to (I5.2p constructed by 
Theorem l5.11 Set z = U\ — u 2 and h = v\ — v 2 with Vj being determined through the divergence-free 
condition with Uj. Then, z and h solve 

rv 

(5.11) dtz + u\d x z + zd x u 2 + v\d y z + hd y u 2 = d 2 z, h = — I d x zdy', 



with z 



y=0 — Z\y=+oo 



0. 



Multiplying the equation for z by e z for some large k, taking integration over T x R + , and 
applying integration by parts, we obtain 

(5.12) f \z\ 2 dxdy + f (k + d x u 2 )\z\ 2 + d y u 2 hz + \d y z\ 2 dxdy = 0. 

2 at JtxK + JtxR 4 
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By the definition of h, we can estimate 

dyUihzdxdy 



J 



TxR 4 



< sup 

t,x 



d y {u2 — U)z(^ J d x zdy'^jdxdy 
y 1/2 d y (u 2 - U)dy^j \\z\\\\d x z\\ 



for some a < 1/2, where || • || denotes the standard L 2 norm on T x R + . Thanks to bounds (|5.4p . 
U2 converges exponentially to U as y — > oo and thus the integral J" R y 1 ^ 2 d y (u2 — U)dy is finite. In 
addition, since the derivatives d x Uj,d y Uj are bounded, by taking k sufficiently large, the identity 
(l5T2j) yields 

d 



(5.13) 



dt 



\z\ dxdy + 



TxR 4 



TxR 4 



1 2 i I |2 

\Z\ + \Zy\ 



dxdy < C\\z x \ 



We will next derive estimates for z y . For this, we take derivative with respect to y to the equation 
for z and multiply the resulting equation by d y z. With noting that z\ y= o = and z yy \ y= o = 
(obtained by setting y = in ([5. lip ), easy computations yield 



1 d 
2~dt 



| z y | dxdy + 



TxR 4 



TxR 4 



-TO I 



'dxdy + 



rl 



TxR 4 



u\d x \z y \ 



+ (Vly + U 2x )\Zy\ 2 + Ul^Z^Zy + U 2 yhyZy + U 2x y Z Zy + f 1 ~ <9j, | | 2 + ^J/^ 



dxdy 



Again, by using the boundedness of Uj x ,Uj xy ,Uj yy , the divergence- free condition h 
similar estimates on the term involving h as above, we easily get 

d 



0. 



-z x , and 



(5.14) 



dt 



TxR 4 



z y \ dxdy < Cl||2;|| + ||zy|| + \\z x 



We note that by using the fact that the derivatives Ujx,Uj xy ,Uj yy are not only bounded, but also 
decay exponentially in y, similar estimates as done above also yield 

d 

(5.15) 



dt 



TxR 4 



y n \z y \ 2 dxdy < C (\\z\\ 2 + \\z y \\ 2 + \\z x \\ 2 



Vn > 0. 



Finally, we may wish to give similar estimates for z x . That is, taking x-derivative to the equation 
for z, testing the resulting equation by z x , and using the boundary condition z x \ y= o = 0, one may 
get 

d 



(5.16) 



dt 



| z x | dxdy + 



TxR 4 



TxR 4 



\z xy \ dxdy 



+ 



TxR 4 



(U\ x + U2 X )|Z X | 2 + U2xx zz x + Vi x Z x Z y + U2yh x Z x + U2 X yhz x 



dxdy 



0. 



However, it is not at all immediate to estimate the term U2 y h x z x in the above identity to yield a 
similar bound as in (|5.14p since h has the same order as z x by its definition (see (|5.1ip ). 

Therefore, we shall derive estimates for z x through the equation (15. 2p and the estimates on w 
obtained in Lemma 15.31 First, we recall that u is defined through w by the relation (see, for 
example, Eq. (4.1.52)]): 

fu(t,x,y)/U(t,x) j 

— n 7\ dr1 '- 

W{t, X, Vj) 
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(5.17) 



u r = u 



U x f U / U W x , 

U J w z 



for it and w being solutions to (jl.3p and ()5.2|) . We apply this expression to ui,wi and U2,W2, 
respectively and derive an estimate for z x = u\ x — u 2x - In regions where u\ > 112, it will appear to 
be convenient to estimate follows: 



Ux\ <C 



(5.18) 



+ \wi 

Whereas in regions where u\ < U2 we estimate 



z\ + \w\{t, X, Ui/U) — W 2 (t, X, U 2 /U)\ 
U2/U 

— T (t,x,T] )dr) +\wi 
ui/U w l 



/•U2/U 


w 2x 






Jo 





{t,x,r}')drf 



u 2 /U 



W\ x W 2x 



UK 



{t,x,rf)drf 



(5.19) 



\z x \ <C 

+ \W2 



z\ + |ioi(i, x, ui/U) - w 2 (t, x, u 2 /U)\ 

u 2 /U 



rui/U 


Wlx 


Jo 


w\ 



(t,x,r]')dr] f 



ui/U 



-¥r(t,x,r)')dr]' 



+ \ W2\ 



Wlx w 2x 



>i'~> 



{t,x,rf)drf 



From the definition Wj(t,x,Uj/U) = d y Uj(t,x,y), we have \wi(t, x, u\/XJ) — w 2 (t,x,U2/U)\ = \z y \ 
Also, note that \wj x /wj\ is uniformly bounded. We have 





w jx 






Jo 


w] 



(t,x,r)')dr)' < C 





1 


Jo 


Wj 



(t,x,Tf)drf = Cy, 



and 



\Wn 



U 2 /U 



ui/U w 



7/1 ■ 

^(t,x,r/)drf <C\ Wj 



U2/U y 

1 7 dr l' 

1 — n 



< c 



Wi 



Wi 



l-m/u i-u 2 /u 



\U\ - U2\ 



lui/U i-lf 

Now, if ux > u 2 , we use the estimate (|5.18p and the fact that \wj\ < C(l — Uj/U). We thus obtain 

\wi\ \wi\ \wi\ 



+ 



< 2- 



< C. 



l-m/U I-U2/ZJ- l-ui/V 

Similarly, if u\ < 112, we use (|5.19p and replace w\ by w 2 in the above inequality, leading to the 
similar uniform bound. This explains our choice of expressions in (|5.18|) - fj5.19j) . By combining 
these estimates, the second and third terms in (|5.18|) when u\ > U2 and in (|5.19|) when u\ < U2 are 
bounded by 

C{\z\+y\z y \). 



Finally, we give estimates for the last term in inequalities (|5.18p and (|5.19p . Using the estimates 
on w,w x , we have 



Wlx W 2a 



w\ 



< C \W\ X - W 2x \ , C \ W l- W 2\ 



[1 - T] 



,'\2 



(l-7/)2 ' 



vVe(o,i), 



^There is an unfortunate typo in [6j Eq. (4.1.53)] where the integral in (|5.17[) was /J* ^-(t,x,r/')dri' 
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which together with the standard Holder inequality implies that 









/ 




/ 


/TxK- 




Jo 



u k /U 



— t )(^ x ^V)dv 



dxdy 



< c 


sup / 




t,X JR + 


< c 


sup / 






< c 


/ 




JTx[0,l] 



dy 



Tx[0,l] 
Tx[0,l] 



\wix - W 2x \ 2 \wi - W 2 \ 2 



{i -if {i -if 

dxdrj' 



dxdrj' 



w\ x ~ w 2x \ 2 \wi - w 2 \ 2 



\wix-W 2x \ 2 , \wi-W2\ 2 



+ 



(I -if {I -if 



{i -if {i -if 

dxdi , 



for some (3 < 3 satisfying (3 — (3)02 < Q\. 
Thus, we have obtained 



(5.20) 



< c 



\4h + \\yz y \\ 2 L 2 + 



Tx[0,l] 



\Wl x - W 2 g 

(l-r,f 



+ 



\Wl — W2\ 



(l-7/)/> 



dxdrj 



for some f3 < 3. Now, applying Lemma 15.31 into (|5.20p . we then have the following estimate: 



(5.21) 



| %X || £2 < C 



kll 2 2 + lly^llia + 



Tx[0,l] 



kla - ^2x| 2 , \lVi - W 2 \ 2 



(i- v y 



+ 



(i- v y 



(0, x, rj)dxdr] 



Combining this with estimates (|5.13p . (|5.14p . and (|5.15p and applying the standard Gronwall's 
inequality, we easily obtain 



(5.22) \\zf Hl (t) < C(T) [\\z \\ 2 m + \\yzoy\\ 2 L 2 + 
where we have denoted zq = uqi — uq2- 



Tx[0,1] 



Wlx ~ W 2x \ 2 \W1 ~ W 2 \ 2 



(l-7]f {I-Vf 



(0, x, i])dxdr) 



Note that ||y2oi/||f,2 < II^^Oj/H 2 ^- It thus remains to express the last estimate in terms of initial 
data uoi and uo2- We note that for r\ = ui(0,x,y)/U(t,x), 

\wi — U>2|(0,X,f7) < \wi(0, X, U\/U) — ^2(0, X, U2/U)\ + \w2{0,X,Ul/U) — W2{0,X,U2/U)\ 

< \d y (ui - u 2 )(0,x,y)\ + I^^Hni - u 2 \{0,x,y). 

In addition, for 77 = «i(0, x, y)/U(t, x), assumptions on initial data (see (O)) gives (1 — rjf 1 < Ce d2y 
and \rj y \ = \d y Uo%/U\ < C(l — uq\/U). Thus, we can make change of variable 77 back to y and 
estimate 



/ 



\Wl — W2\ 



Tx[0,l] iX-^f 



(0,x,ri)dxdn < C e ( P- 1 ^ y (\d y (u 01 -u 2)\ 2 + \u 01 -u 2\ 2 ){x,y)dxdy 



Txlft, 



< c\\e^-^y/ 2 ( U01 - U0 2)\\ 2 H1 . 

Similarly, we have 

\wi x -w 2x \(0,x,ri) < \d x d y (u i - u 02 )\(x,y) + C\d x (u i - u 02 )\(x,y), 

and thus 

/ ]Wl ( x - W2 f (0,x, V )dxd V < C\\e^-^y/ 2 (u 01 -u 2)\\ 2 H2 . 

Ax [0,1] ^-"nf 
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Putting these into (|5.22|) . we have obtained 

(5.23) Il(«i-«2)(t)[&i < C\\e a v(uox-u 02 )\\ 2 H2 , 

for a = {f3- l)9 2 /2. Theorem [___ thus follows. □ 
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